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Abstract 

An algebraic model based on Lie-algebraic and discrete symme- 
try techniques is applied to the analysis of thermodynamic vibrational 
properties of molecules. The local anharmonic effects are described by 
a Morse-like potential and the corresponding anharmonic bosons are 
associated with the SU(2) algebra. A vibrational high-temperature 
partition function and the related thermodynamic functions are de- 
rived and studied in terms of the parameters of the model. The idea of 
a critical temperature is introduced in relation with the specific heat. 
A physical interpretation of a quantum deformation associated with 
the model is given. 

1 Introduction 

The algebraic approach has been used successfully in molecular physics and 
has led to new insights into the nature of complex many body systems 
jl], |||. In the framework of the algebraic method, the Hamiltonian of a 
given system is written as an algebraic operator using the generators of the 
appropriate Lie algebra. All other operations in the model are algebraic op- 
erators, unlike the differential operators in standard wave mechanics. The 
technical advantage of the algebraic approach is the comparative ease of al- 
gebraic operations. Equally important, however, is the conclusion derived 
from comparison with experiment, that there are generic forms of symmetry- 
adapted algebraic Hamiltonians and that entire classes of molecules can be 
described by generic Hamiltonians where the parameters vary in a systematic 
fashion for different molecules. In its initial stage of development 0]- p3] , the 
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algebraic approach has sought to show why and how it provides a framework 
for the understanding of large- amplitude anharmonic motion. The anhar- 
monicities are introduced by means of dynamical groups that correspond to 
anharmonic potentials and which constrain the total number of levels to a 
finite value. Later on, the SU(2) models [15|- p(J combined Lie algebraic 



techniques, describing the interatomic interactions, with discrete symmetry 
techniques associated with the local symmetry of the molecules. Recently, 
a clear-cut connection could be established between the Moise-SU(2) ap- 
proach and the traditional potential energy surface methods P, p), 0]. 

The algebraic model [|], exploits the isomorphism between the SU(2) 
algebra and the one-dimensional Morse oscillator, 

where D is the depth of the potential well, d is its width, x is the dis- 
placement from the equilibrium and \i is the mass of the oscillator. The 
one-dimensional Morse Hamiltonian can be written in terms of the genera- 
tors of SU(2), 

H M = j(n 2 - 4Jf) = yGM- + J-J+- N), (2) 

where A is a constant dependent on the parameters of the Morse poten- 
tial. The eigenstates, ||_iVj , ^) , correspond to the U{2) D SU(2) symmetry- 
adapted basis, where N is the total number of bosons fixed by the potential 
shape, and v is the number of quanta in the oscillator, v = 1,2,... , [yj . 
The value of TV" is dependent on the depth D and the width d of the Morse 
potential well |, §, [Of, 

The parameters A and iV are related to the usual harmonic and anharmonic 
constants u) e and x e u e used in spectroscopy ||, 

/ 2D \ 2 

Ue = A(N + l)=h\J£) , (4) 

_ , n 2 
XeUJe ~ Wd- 
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The anharmonic effects can be described by anharmonic boson operators 



J+ f t J- . N t 

= , b> = —= , v = J z 

N VN 2 2 



(5) 



where v is the Morse phonon operator with an eigenvalue v. The operators 
satisfy the commutation relations, 



b, v 



b ,v 



b,b 



1 



2v 
~N 



b, V 



(6) 
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The harmonic limit is obtained when N — ► oo, in which case 
giving the usual boson commutation relations. 

The Morse Hamiltonian can be written in terms of the operators b and 



H M ~~(bP+Pb 
which corresponds to vibrational energies 



1 v J 

: , - huj [ v + - - — ) , v = 1,2, ... , 



N 
~2 



(7) 



(8) 



where ujq is the harmonic oscillator frequency. Thus, the spectrum of the 
Morse potential leads to a deformation of the harmonic oscillator algebra. A 
more detailed relationship between the Morse coordinates and momenta and 
the SU(2) generators can be derived through a comparison of their matrix 
elements Jl2]| and through the derivation of raising and lowering operators for 
the Morse potential | |13| |. Note that for an infinite potential depth, N ^oo, 
the Morse potential cannot be distinguished from the harmonic potential. 
To sumarize, the algebraic anharmonic model has been developed to analyze 
molecular vibrational spectra |2||-pC|]. It provides a systematic procedure for 
studying vibrational excitations in a simple form by describing the stretching 
and bending modes in a unified scheme based on SU{2) algebras which 
incorporate the anharmonicity at the local level. 

The deformation of the harmonic oscillator algebra, associated with the 
Morse potential, has also been derived using a quantum analogue for the 
anharmonic oscillator 
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We have described the anharmonic vibrations 
as anharmonic g-bosons using a first-order expansion of the quantum de- 
formation. We have thus proposed a physical interpretation of quantum 
deformation in the framework of the algebraic model. 
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The aim of this paper is to apply the algebraic approach to the vi- 
brational high-temperature thermodynamics of diatomic molecules and to 
obtain the basic thermodynamic functions in terms of the parameters of the 
algebraic model. This paper can be considered as a first step in the direction 
of incorporating anharmonicity and the finite number of bound vibrational 
states into the thermodynamic description of molecular systems. The rest 
of the paper is organized as follows. In Section 2, we derive a Morse-like 
vibrational partition function for high temperature and study its properties. 
In Section 3, the vibrational partition function is used to derive the basic 
thermodynamic functions, such as the mean vibrational energy, specific heat 
and free energy. The idea of critical temperature is introduced in relation 
with the specific heat. In Section 4, the mean number of anharmonic bosons 
is obtained. The concept of maximal temperature of the anharmonic vibra- 
tions is discussed. The g-bosonic deformation of first order is considered. 
It is shown that this quantum deformation is related to the shape of the 
anharmonic potential well and the fixed number of anharmonic bosons. The 
results are applied to the diatomic molecule H 1 C1 35 . These results must be 
combined with the translational and rotational thermodynamic functions in 
order to compare with experiment, as discussed in reference J21[ |. 



2 Vibrational Partition Function 

The vibrational partition function of a diatomic anharmonic molecule is 

[N/2] 

Zn=J2 e ~ P£v ( 9 ) 

v=0 

where (3 = 1/ksT, the vibrational energies e v are given by equation (g) 
and N is the fixed total number of anharmonic bosons. Introducing new 
parameters, a = , iVo = [^J and I = Nq — v, the exact value of 
vibrational partition function can be written as, 

N 

Z N = e- a J2 e ~ fr ^ Ni ~ l2 " > - ( 10 ) 

1=0 

At high temperatures T, for Nq large and a < 1, the sum can be replaced 
by the integral, 



/aNo 

— at Nn+1) 
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Figure 1: Vibrational partition function Z^q as a function of a. The double- 
dashed line represents the exact representation. 



where s = yj^fcl- This integral can be evaluated exactly in terms of the 
error function, erf iU/olNq) (as defined in (23|]), 



Z N = \^ ^e-^%rU(^I%) . (12) 

Equation (^) represents the high-temperature value of the vibrational 
partition function in the Morse-like spectrum [24, ^5], 26]. The partition 
function is expressed in terms of the parameters of the algebraic model Nq 
and a. The dependence on the temperature is given by a, 

a = — ; — ^ = (13) 
2k B T 2T y ' 
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where = hioo/kB is the usual characteristic vibrational temperature of 
the molecule. The contributions of the anharmonic vibrations are essential 
in the high-temperature region for T > 0, where T = corresponds to 
a = 0.5. 

When N$ — > oo, the harmonic limit of the model is obtained, 

Noe~ a e~ a T e 

Z^, ~ ~ == — e~2T (14) 

00 2aiV -l 2a 1 J 

which coincides with the harmonic vibrational partition function of a di- 
atomic molecule at high temperatures. The expression for the partition 
function (^) can be generalised to polyatomic molecules, by combining the 
present results with the use of a local-mode model where each interatomic 



potential is of the Morse form [14] 



The diatomic molecule H 1 C1 35 is considered as an example. Using the 



values of x e and uj e [21] for the zero lines of this molecule and equations 
(||) and (||) we obtain the total number of anharmonic bosons, fixed by the 
shape of the Morse potential, N = 56, and the total number of quanta in the 
oscillator, Nq = 28. The depth of the Morse potential is D = 5.32eV, and 
the width is d = 0.57 x 10~ 10 m. The characteristic vibrational temperature 
of the molecule is = 4300K 



Substituting the value of iVo = 28 in equation (12) we can calculate 
the partition function, Z56, for the molecule H 1 C1 35 as a function of the 
parameter a. The graph on Figure |l] represents the partition function Z§q 
given by equation (|l2|) for the values of the parameter a between and 1 
(solid line). The exact partition function from equation ( |l0|) is given for 
comparison (double-dashed line). It is clear that in the region < a < 
1, the integral approximation is in a very good agreement with the exact 
representation and does not change the value and appearance of the partition 
function. The comparison between the function Z56 (solid line) and the 
harmonic limit Z^ (dashed line) is given on Figure ^. The finiteness of Z§§ 
in the high T limit is linked, of course, with the finite number of states in the 
Morse potential. A more realistic description in the high T region requires 



the introduction of the continuum states of the Morse potential [28]. 

Having the partition function Zn in terms of the parameters of the 
algebraic model, we are now in position to derive the corresponding ther- 
modynamic functions. An algebraic approach has been used in pj| to study 
the thermodynamic properties of molecules. However, the partition func- 
tion in |29| uses an approximation of the classical density of states, while we 
have derived an explicit function in terms of the parameters of the algebraic 
model. 
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Figure 2: Vibrational partition function Z§§ (solid line) and the harmonic 
limit Zqo (dashed line) as functions of a. 



3 Thermodynamic Vibrational Functions 

3.1 Mean Vibrational Energy 

The mean vibrational energy is given by 

d huj dZ N 

Taking into account that 



dZ N _ Z N , N e-« 

-^--^- {No + 1)Zn + ^t> 



(15) 



(16) 
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Figure 3: Mean vibrational energy U^/frwo as a function of a. The dashed 
line represents the harmonic limit Uoo/Tuvq. 



we obtain the following expression for the mean vibrational energy in terms 
of the partition function Z^, 

2 V 2a 2aZ N J 

Substituting Zn with equation (^) gives the following expression for the 
mean energy, Un, in terms of the parameter a, 

2 2a V «7T erf % (VaiVo) / 



The harmonic limit is obtained from equation (|17D, when Nq — > oo and Zn 
is given by (PI), 

C/oo ~ -^-(1 + -) = — + A; S T. (19) 
2 a 2 

This is the classical mean energy of a diatomic molecule at high tempera- 
tures. 

The graph on Figure ^ represents the mean vibrational energy, U^q / ftwo , 
of the molecule H 1 C1 35 for values of a, < a < 1. The high-temperature 
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region corresponds to a < 0.5. For comparison the graph of the harmonic 
limit Uoq (dashed line) for No — > oo is also given. As already mentioned in 
the previous section, the finiteness of Un in the high T limit is a result of 
the finite number of states in the Morse potential. 



3.2 Specific Heat 

The vibrational part of the specific heat is, 

dU N huj dU N 

C N = — = -^-^ — . (20) 

Substituting Un with equation ( |i~7| ) and using equation (|T6|), we obtain 

This equation represents the vibrational specific heat in the algebraic model 
in terms of the partition function Zjy. Substituting Zn in equation (|2~l|)with 
the expression fll2|) , we obtain the dependence of the specific heat Cat on 
the parameter a, 



2 f f erf* (V^No~) \ 2 V 7T erf % (VaiVo) J 

(22) 

It is clear from the relation ( p2|) that all anharmonic contributions to the 
vibrational part of the specific heat depend on the parameter a and hence on 
the temperature T. When A^o — > oo, the harmonic limit of the model gives 
the vibrational specific heat of a diatomic molecule at high temperatures, 

Coo ~ ks- (23) 

Figure || represents the dependence of the vibrational specific heat, C^q /ks, 
on the parameter a, < a < 1, for the molecule E^Cl 35 (solid line). The 
graph of the harmonic vibrational specific heat of a diatomic molecule (see 
.g. |27j), Charm/ks, is also given (double-dashed line) as well as the 
harmonic limit, Coo/fog (dashed line), where 

Charm = k B ^C? — -5. (24) 

(e Aa — 1 
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Figure 4: Vibrational specific heat C^/ks (solid line) as a function of a. 
For comparison, Charm/ ks (double-dashed line) and C^/ks (dashed line) 
are also given. 

The effects of the anharmonicity are strongest for values of a < 0.5, where 
a = 0.5 corresponds to the characteristic vibrational temperature © (6 = 
4300K for H^Cl 35 ). 

The graph shows an anomaly in the dependence of the vibrational specific 
heat on the parameter a (temperature T). The specific heat has a maximum 
for a value of a = olq which corresponds to a temperature, T = Tc- We 
shall call this temperature the critical temperature and the corresponding 
parameter ac the critical parameter. The anomaly of the specific heat 
is again a result of the finite number of states in the system. The specific 
heat increases with the temperature as more anharmonic bosons are excited. 
The maximum is reached when the latter occupy the energy state with 
v = Nq = 28. (Note that the shape of the curve is similar to the shape 
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associated with the Schottky anomaly of the specific heat of a two-level 
system p7|| ). 

The maximal vibrational energy is given by the equation (||) by substi- 
tuting v with its maximal value Nq. Thus, 



while the minimum energy is, 



( ^ • * (2o) 



e = Ijhwo . (26) 



Thus, 



Ae = e max - eo = ^hu) N (27) 



This gives Ae = = 5.1877eV for the molecule H 1 C1 35 . Comparing 

Ae with the dissociation energy of the molecule DE = 4.4703eV |3(J, we 
can conclude that at the temperature T = Tc, Ae > DE and some of the 
molecules might have started to dissociate while others may still be in stable 
molecular states. Our model, in its present form, does not account for the 
effects of the dissociation. In addition, this simple version of the model does 
not yet include the contributions of the translational and rotational degrees 
of freedom which at temperatures close to Tc may be substantial. The 
critical temperature Tc can be considered as a temperature above which 
the model is no longer valid in its current form and other effects take place, 



.g. dissociation [28]. 



We have studied the behaviour of the specific heat with respect to the 
combined parameter o.Nq. The graph of Cjv /&b as a function of o.Nq shows 
a similar anomaly (Fig. |5|). 

Solving numerically the equation ^n ) = ^ with respect to the com- 
bined parameter aNo, we have found a root, oicNq = 6.1332642. Thus, 
the critical value ac decreases as the number of fixed anharmonic bosons 
increases, 

6.1332642 . . 

a c = (28) 

No 

When iV"o — > do, ac — * and the anomaly of the specific heat disappears, 
which is in agreement with the harmonic limit of the model. 

For the molecule I^C 35 , N = 28 which gives a value for ac = 0.219. 
Thus, the value of the critical temperature for this molecule is Tc = 9815-ftT. 
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Figure 5: Vibrational specific heat CV/fce as a function of ctNa 



3.3 Free Energy 

The free vibrational energy in terms of the partition function Z^ris given by 



F N = --\nZ N 





(29) 



Substituting with equation (12) gives the free vibrational energy in the 
algebraic model at high temperatures, 



F N 



hu>o 



I m2 + J_ln f-£L.\ + (N + 1) - -In (erf i (aN )) 
a 2a \ tti\q J a 



(30) 



When iVo — * oo, using expression (|1J) in equation fl29|), we obtain the 
classical harmonic result for the free vibrational energy at very high tem- 
peratures, 



Foo ~ k B T\n2. 



(31) 



4 Anharmonic Bosons 
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4.1 Mean Number of Anharmonic Bosons 

The mean vibrational energy in the anharmonic model can be written in 
terms of mean number {un) of anharmonic quanta, each with energy fuvo, 

U N = hu ({v N ) + \) (32) 



Substituting Un by equation (|i7|), we obtain (vn) in terms of the partition 
function Zjy, 

No 1 N e- a 

M = T + Ta ~ W (33) 

Using the expression (^) in equation (^), we obtain the high-temperature 
value, 



N 1 / No e aN ° 
{UN) = ^ + ^-\l4^ern(V^No.) ^ 

The harmonic limit is obtained from equation ( pip when iVo — > oo and 
Zjy is given by the expression (14), 



The graph of the function (u^q) for the molecule H^C 35 is given on Figure 
||. The dashed line represents the harmonic limit (z^oo)- The same reasons 
apply to the finiteness of (v^e) as those discussed for the partition function 
in Section 2. 



4.2 Maximal Temperature 

The maximal vibrational energy (^) is obtained when v = No. The maximal 
mean vibrational energy is given by equation (^) when vn = (^v)maxj 

Umax = ftwo(W)max + ~) • (36) 

Comparing the equations ( |25|) and ( |36|) gives the maximal mean number of 
anharmonic bosons, 

(fjv)max = "TT ( 37 ) 
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Figure 6: Mean number of anharmonic bosons v^q as a function of a. 



(which according to the graph on Figure || is for a value a — > 0) . Substituting 
(i'jv) by equation ( |34] ) and simplifying gives, 

2 [a^_ e aN = erf ■ (3g) 

The numerical solution of the above equation has a root aiVo— »0. As iVo is 
a fixed number, this solution leads to a — > 0. In terms of the temperature, 
we obtain T max — ► oo. The result shows that in practice the system does 
not reach a maximal temperature, which shows the need of incorporating 
the continuum into the partition function (||). 



4.3 Quantum Anharmonic Bosons 

In p2l, we have shown that the anharmonic bosons b, 



can be obtained as an approximation of g-bosons [31, |3 




from Section 1 
3|. The g-bosons 
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are defined by the following commutation relations: 



a, af] = q n , [n,a] = —a , [n,a)\ = a) (39) 



where the deformation parameter q is in general a complex number |32|| . As 
shown in [22], the anharmonic commutation relations (||) can be recovered 



for real values of the deformation q close to 1, q < 1, and a linear expansion 
of q in terms of a parameter p, p = 1/(1 — q), 

Tl 

q n = l-- . (40) 
p 

If we now substitute the aproximation for q n from equation (|4(]) in the 
commutation relations ( |39[ ) and identify the parameter p with N/2, h with 
v and the creation and annihilation operators a, a) , with 6, &t respectively, 
we recover the SU{2) anharmonic relations @. 

The form (||) of the SU (2) commutation relations can be considered as a 
deformation of the usual (harmonic oscillator) commutation relations, with 
a deformation parameter p = N/2. This gives a possible physical realization 



for the quantum deformation obtained in [22]: the quantum deformation pa- 
rameter p is the fixed number Nq of the anharmonic bosons in the oscillator. 
Using the relation between the fixed number of anharmonic bosons N and 
the characteristics of the Morse potential (||) we arrive to the conclusion that 
the quantum deformation is also determined by the depth, the width and 
in general the shape of the Morse potential well. For the molecule H 1 C 35 , 
p = 28 which gives q = 27/28. 

Now, substituting iVo = p in the expressions for the partition function 



(12), mean energy (|Tq) , specific heat (f2j), free energy (p0|) and mean num- 
ber of anharmonic bosons (|33|), we obtain the thermodynamic properties of 
diatomic molecules in terms of the deformation parameter p. Equation ( |28| ) 
gives the relation between the quantum deformation parameter and the crit- 
ical parameter ac (critical temperature Tq). For large values of p, (q— 
the classic harmonic case is restored. 



5 Conclusion 

We have studied the vibrational thermodynamic properties of diatomic molecules 
which at high temperature may strongly depend on its anharmonicity. We 
have derived and studied the vibrational partition function and the related 
thermodynamic functions, such as mean vibrational energy, specific heat 
and mean number of anharmonic bosons, in terms of the parameters of the 
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model. We have also shown that it is possible to interpret these results 
in terms of a quantum deformation, related to the shape of the Morse po- 
tential, and which is associated with the fixed total number of anharmonic 
bosons, so that the thermodynamic properties of diatomic molecules depend 
on the corresponding quantum deformation parameter. We believe that this 
paper constitutes a first step in the description of thermodynamical prop- 
erties of diatomic molecules which in principle can be simply generalized 
to polyatomic molecules. We are currently studying the introduction of 
the continuum into the description, in order to take into consideration the 



transition to dissociation 28 
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